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1. Introduction 

The three-dimensional incompressible Euler equation 

— + (u ■ V)u = —VP , divu = 
at 

describes the motion of an inviscid incompressible fluid. A steady solution of the 
Euler equation is a time-independent vector field u that satisfies the above equation 
for some pressure function P. Equivalcntly, the vector field u satisfies 

u A curl u = V-B , 

where B := P+\u\ 2 /2 is the Bernoulli function. In fluid mechanics, the trajectories 
of the velocity field u and of its curl (i.e., the vorticity) are called stream lines and 
vortex lines, respectively. 

Since Lord Kelvin, steady solutions of the Euler equation having knotted and 
linked stream (or vortex) lines have attracted considerable attention, partly be- 
cause of their connections with turbulence and hydrodynamic instability [3J 112j . 
Because of the connection of the steady Euler equation with variations of the en- 
ergy functional under volume-preserving diffcomorphisms, the presence of periodic 
trajectories of prescribed link type has also been considered to derive lower bounds 
for the energy functional in certain natural classes of vector fields. These bounds 
can actually be expressed in terms of the helicity and asymptotic linking number 
of a fixed vector field of the class (cf. e.g. [3 [25] and references therein) . 

Despite the large body of related literature, there are surprisingly few results of 
existence of steady solutions with linked stream lines, and in fact our knowledge 
is mostly based on numerical simulations and exact solutions. In the mid eight- 
ies, Moffatt [20] proposed an appealing heuristic approach to the construction of 
steady solutions with stream lines diffcomorphic to a given link. Unfortunately, 
this approach requires global well-posedness for a magnetohydrodynamic system of 
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PDEs including a Navier-Stokes equation (which, as is well known, is very far from 
obvious) and makes strong topological assumptions on the asymptotic behavior of 
the solutions as t —> oo. 

To our best knowledge, there are only two existence results of steady solutions 
with prescribed topological properties. Firstly, using contact topology, Etnyre and 
Christ (i showed that, given a certain link, there exists a (possibly incomplete) 
Riemannian metric in M 3 adapted to it such that the Euler equation in this metric 
admits a steady solution having this link as a set of periodic stream lines. Consid- 
ering the Euler equation for arbitrary Riemannian metrics is key in order to make 
the problem amenable to a purely topological approach, so this method cannot be 
used to derive existence results for a fixed (e.g., Euclidean) metric. 

The second result, by Laurence and Stredulinsky |16j . utilizes variational tech- 
niques to analyze steady solutions of the Euler equation with axial symmetry. This 
symmetry makes the problem effectively two-dimensional and allows to formulate 
it in terms of the scalar stream function. In this case, the authors prove the ex- 
istence of continuous, energy-class solutions in a bounded domain with level sets 
of a given topological type. This approach cannot be modified to prescribe stream 
lines because it provides some (weak) control over the poloidal component of the 
velocity field but does not yield any information on its axial component. Moreover, 
at most the so-called torus knots could be possibly realized as stream lines of an 
axisymmetric solution. 

In this paper we aim to fill this gap by proving the existence of steady solutions 
in R 3 with stream or vortex lines of prescribed link type. Here by a link we mean 
a disjoint union of tame knots, that is, of smoothly embedded circles in R 3 . The 
steady solutions we construct belong to the class of (strong) Beltrami fields (TTJ 
Section 2.3], i.e., they satisfy 

curl u = Xu 

with A being a nonzero constant. Beltrami fields also play a preponderant role in 
magnetohydrodynamics, where they are called force-free fields. It should be noticed 
that the stream lines of a Beltrami field are also vortex lines. 

The motivation to consider Beltrami fields to produce steady solutions with 
stream or vortex lines of arbitrary link type comes from Arnold's and Etnyre and 
Christ's studies of solutions to the steady incompressible Euler equation when the 
Bernoulli function B is analytic and nonconstant [TJ [5] . Under certain dynamical 
assumptions, these authors showed that in this case the only knots that could 
possibly be stream or vortex lines of a steady solution to the Euler equation in a 
bounded domain are those of Wada type; this result also holds in the whole R 3 
provided that the solution vector field satisfies the estimate \u(x) \ + | curlu(x)| ^ 
A + B\x\ for some constants A and B. That a similar result does not hold for 
Beltrami fields was well known through the analysis of particular instances of ABC 
flows [5] . Our main result in this paper is that actually there are no restrictions on 
the knot types that can be realized as stream lines of a Beltrami field: 

Theorem 1.1. Let Lcl 3 be a possibly unbounded, locally finite link. Then for 
any real constant A ^ one can transform L by a C°° diffeomorphism <f> of R 3 
arbitrarily close to the identity in any C r norm, so that $>{L) is a set of stream 
lines of a Beltrami field u, which satisfies curlu = Xu in R 3 . 
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The proof of this theorem is given in Section [3j the demonstrations of several 
intermediate results being postponed to Sections [4|j7| An informal guide to the 
proof of the theorem is presented in Section [2] The proof combines techniques from 
ordinary and partial differential equations with differential topology. While it is 
not easy to find parts of the proof purely relying on methods from just one of these 
areas, the basic philosophy is to use differential topology and ordinary differential 
equations (as opposed to Laurence Strcdulinsky) to gain control on the trajectories 
of various auxiliary vector fields, and partial differential equations to relate these 
auxiliary vector fields to the Beltrami equation (as opposed to Etnyre-Ghrist). 

Theorem 1 1 . 1 1 provides a rigorous proof of Moffatt's heuristic insight [2D] that any 
knot should be realizable as a stream line (or vortex line) of a steady solution of the 
Euler equation. It also furnishes a positive answer to Etnyre and Ghrist's question 
of whether there exists a steady solution in R 3 having linked stream lines of all 
isotopy types; details and further applications will be given in Section |U 

It should be noticed that any vector field satisfying curlu = Xu has infinite 
energy, since in this case Am = — X 2 u and the Laplacian does not have any L 2 
eigenfunctions in ]R 3 . While there are L°° solutions of the equation curl u = Xu 
in M 3 , such as the ABC flows T7\ . the solutions given by Theorem are not 
granted to be bounded either. Incidentally, the related problem of whether there are 
any nonzero L°° steady solutions to the three-dimensional Navier-Stokes equation 
is wide open; the fact that no bounded solutions exist either in two dimensions or 
in the axisymmetric case has been established recently [13] . 

2. Strategy of proof 

Leaving technicalities aside for the moment, the strategy of the proof of the 
main theorem is the following. We take a connected component L a of the link 
L, which is a smooth knot, and perturb it a little if necessary to make it real 
analytic. We then embed the perturbed knot into an analytic two-dimensional 
strip and construct a vector field on this strip having the latter knot as a stable 
hyperbolic limit cycle. Next we extend this vector field on the strip to a local 
Beltrami field v a , defined in a tubular neighborhood of the cycle. By construction, 
the field v a has a periodic orbit diffeomorphic to the link component L ai which is 
granted to be hyperbolic by the choice of the field on the strip and the fact that 
the local Beltrami field v a is divergence-free. The above procedure is then repeated 
for all the components of the link. To conclude, we show that there exists a global 
Beltrami field u that approximates each local field v a ; the field u will then have a 
stream line diffeomorphic to each link component L a because of the robustness of 
the limit cycle of v a . 

Three main difficulties arise when trying to carry out this program. In the 
first place, it seems natural to define the local Beltrami field v a as the solution to 
some Cauchy problem for the curl operator, but the fact that the symbol of curl 
is an antisymmetric matrix (and thus there are no non-characteristic surfaces of 
curl) precludes a straightforward application of the Cauchy-Kowalewski theorem. 
In the second place, the approximation result connecting the global Beltrami field 
u and the local fields v a should be sufficiently fine to effectively deal with the 
noncompactness of the link L. Thirdly and related, the vector field on the strip 
used as "Cauchy data" to define the local Beltrami field v a must be delicately chosen 
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in order to ensure that the limit cycle of v a diffcomorphic to the link component 
L a is "robust enough" to be preserved by the latter approximation. 

The proof of Theorem which we present in the following section, is divided 
into five steps. Let us very briefly comment on the purpose of each step: 

(i) First of all, in Step 1 we present an abstract local existence theorem (The- 
orem 3.1) for the equation curlw = Xv with data v\s = w on a surface E 
that will be subsequently applied to define the local Beltrami field v a in 
terms of a prescribed vector field on a strip at the link component L a . To 
circumvent the lack of ellipticity of curl (or of its counterpart *d acting on 
1-forms, with * the Hodge operator), we resort to an indirect argument 
with the Dirac-type operator d+d* acting on the space of differential forms 
of mixed degrees. Here d* is the codifferential; the operator d + d* is then 
elliptic (thus allowing us to invoke the Cauchy-Kowalewski theorem), but 
it does not preserve the space of 1-forms. Ultimately, this reflects itself in 
the necessity of imposing some technical conditions on the vector field on 
the strip that will serve as Cauchy data. 

(ii) Appropriate Cauchy data (that is, the vector field on the strip at each link 
component L a previously mentioned) are carefully constructed in Step 2 



(Lemma 3.4). Each field is required to satisfy the technical conditions 
demanded by the existence theorem proved in Step 1 and, besides, to have 
a stable hyperbolic limit cycle diffeomorphic to the link component L a . 
This construction is considerably simplified by the introduction of a local 
coordinate system adapted to the link component, 
(iii) The choice of Cauchy data on the strip we made in Step 2 ensures that the 
associated local Beltrami field v a obtained using the existence theorem 
of Step 1 also has a periodic orbit diffcomorphic to the link component 
L a . In Step 3 we show that this periodic orbit is robust under suitably 
small perturbations of the vector field because our choice of Cauchy data 
guarantees that there is a stable manifold tangent to the strip and, the 
field v a being divergence-free, there must also be an unstable manifold, 



thereby granting that the periodic orbit is hyperbolic (Theorem 3.5) 



(iv) In Step 4 we prove an abstract result on approximation of local Beltrami 



fields by global Beltrami fields (Theorem 3.6 ) adapted to the kind of per- 
turbations that were allowed in Step 3. This result relies on a better-than- 
uniform approximation theorem for the auxiliary scalar elliptic equation 
Af = —A 2 /, which in turn hinges on an iterative scheme that employs the 
Lax-Malgrange theorem, 
(v) To conclude, in Step 5 we use the latter approximation theorem to con- 
struct a global Beltrami field u that approximates each local Beltrami field 
v a in a neighborhood of the link component L a and show that therefore 
the field u has a set of periodic trajectories diffeomorphic to the link L on 
account of the hyperbolic permanence theorem. 



3. Proof of the main theorem 

In this section we shall present the proof of Theorem |1.1| which, as we have just 
mentioned, is divided into five steps. We will henceforth assume that the parameter 
A that appears in the Beltrami equation curlu = \u is positive (which can be 
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accomplished by mapping x i— > —x if necessary) and that all the diffcomorphisms 
are of class C°° . 

Step 1 (Local existence for the Cauchy problem). A natural way of obtaining local 
solutions of the Beltrami equation with partially prescribed behavior is to consider 
the Cauchy problem 

curlw = Ai>, u|s=u», 
where E is an embedded oriented surface in K 3 of class C w (where, as customary, 
C u stands for real analytic). Therefore, our goal in this step is to establish the 
following abstract existence theorem for the above Cauchy problem, which grants 
the existence of a solution under suitable hypotheses on the vector field w used as 
Cauchy data: 

Theorem 3.1. Let E be an embedded oriented analytic surface in M 3 . Let w be a 
C u vector field tangent to E and let us denote by 7 its associated 1-form. If the 
pullback of the 1-form 7 to the surface E is closed (i.e., dj^^) — with js : E — > R 3 
the inclusion map), then the equation curlw = Xv with Cauchy data v\s = w has a 
unique solution in a neighborhood of the surface E. This solution is analytic. 

The proof of this result is given in Section |4j To understand the content of 
the theorem, it is worth discussing at this point the theorem's hypothesis that the 
pulled-back 1-form ^(7) is closed. To this end, let us take C u local coordinates 
(p^ 1 ^ 2 ) m a neighborhood of the surface E, where p denotes the signed distance 
to this surface [14]. In these coordinates, the Euclidean metric reads 

(3.1) ds 2 = dp 2 + h lJ (p,Oded? , 

and obviously the surface E is the zero set of the coordinate p. We will denote by 
h lJ the inverse matrix of hij and by \h\ :— det(/ijj) its determinant. Let us write 
the vector field v in these coordinates as 

Q Q 

setting a, := hijb 3 , its associated 1-form can be written as 

P^xiPiOdp + a^OdC- 

Let us suppose that the field v satisfies the equation curlw = Xv and is tangent 
to the surface E, as claimed in the theorem. Since the coordinate vector field d/dp 
is orthogonal to d/dt? by the coordinate expression of the metric (cf. Eq. ( |3.1[ )), 
the condition that the field v be tangent to the surface is obviously tantamount to 
demanding that the function y appearing in the coordinate expression of v satisfies 
X(0,0 = for all £. 

It is slightly more convenient to attack the equation curl v = Xv using differential 
forms, in terms of which it can be rewritten as = Xf3, where /3 is the 1-form 
dual to the field v, * denotes the Hodge operator and d is the exterior derivative. 
The action of *d on the 1-form f3 is given in coordinates by 

+i<» i,v '(f-fK 
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Accordingly, the equation *d(3 = A/3 and the condition x(0, £) = imply that the 
coefficients aj(p, £) of the 1-form /3 satisfy 

This means that the pullback j'^(/3) to the surface of the 1-form /3 must be closed. 
Hence the content of Theorem 3.1 is that this condition (which we state in terms 
of the Cauchy data) is not only necessary, as the calculations above show, but also 
sufficient. 



Step 2 (Construction of appropriate Cauchy data). In view of the local existence 
theorem proved in Step 1, our goal in this step is to construct vector fields tangent to 
certain analytic surfaces that satisfy the hypotheses of the aforementioned theorem 
and, in addition, possess certain convenient dynamical properties on these surfaces 
that will be discussed later. Since we will need to work in the analytic category, 
as a preliminary technical result we start by recording that there is a "small" 
diffeomorphism <t of R 3 which transforms the link L into an analytically embedded 
submanifold Q(L). More precisely, 

Lemma 3.2. Given the locally finite link L, we can transform it by a diffeomor- 
phism $ of M. arbitrarily close to the identity in any C r norm, so that Q(L) is a 
real analytic link. 

We omit the proof of this result, which is a straightforward consequence of 
Thorn's isotopy theorem (cf. e.g. [4j Theorem 14.1.1]). Replacing the link L by 
its transformed image <&(£) if necessary, by the above lemma we can (and will) 
assume that the link L is analytic without loss of generality. We will denote the 
connected components of the link by L ai where the labeling index a ranges over an 
at most countable set A. 

In order to apply the Existence Theorem |3.1| later on, we need to introduce an 
appropriate surface S a containing the link component L a . As the normal bundle of 
L a is trivial |19j . a technically convenient way of doing this is via a C w trivializing 
map a whose 0-fiber is L a : 

Definition 3.3. We say that an embedded analytic surface S a is a strip at the link 
component L a if there is a tubular neighborhood N a of the link component and an 
associated C u trivializing map 9 a : iV a — » M? whose 0-fiber is L a and such that 
the surface S Q is given by 8~ 1 ((— 1, 1) x {0}). Obviously the link component L a is 
contained in the strip S a . 

Clearly the surface E a is diffeomorphic to the cylinder S 1 x R, so it will be 
assumed to be oriented. The following lemma provides a way of choosing a vector 
field w a tangent to a strip at a link component L a so that, on the one hand, 
satisfies the hypotheses of the Existence Theorem |3 . 1 1 that we stated in Step 1 and, 
on the other hand, its pullback to the surface has the link component L a as a stable 
hyperbolic limit cycle. This simply means that there exists a neighborhood of the 
link component in the strip E a whose w-limit set along the local flow of the pulled 
back vector field is this component and that the associated monodromy matrix does 
not have any nontrivial eigenvalues of modulus 1. 
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Lemma 3.4. Given a link component L a , one can take a tubular neighborhood 
N a of L a and an associated strip E a such that there is an analytic vector field w a 
in N a tangent to the strip and whose pullback to the strip (i.e., i„(w ) with j a : 
S a — > N a the inclusion map) has the link component L a as a stable hyperbolic limit 
cycle. Furthermore, one can assume that the pullback to the strip E a of the I- form 
associated to the field w a is closed and that tubular neighborhoods corresponding to 
distinct link components have disjoint closures. 

The proof of Lemma [374] is presented in Section [5] One should note that the fact 
that the tubular neighborhoods corresponding to distinct link components have 
disjoint closures will be used in Step 5 to construct the diffcomorphism $ which 
transforms the link into a set of stream lines. 

Step 3 (Stability of periodic stream lines). Given a link component L a , we showed 



in Step 2 (Lemma 3.4) that one can take an appropriate surface E a containing L a , 
which we called a strip, and a vector field w a tangent to the strip which has certain 
convenient dynamical properties and satisfies the requirements of the Existence 
Theorem |3. 1 1 proved in Step 1. Therefore, there is a unique analytic local Beltrami 
field v a , which satisfies the equation curl?; a = Xv a in a small enough tubular neigh- 
borhood N a of the link component, whose restriction to the strip is w a |s a = w a . 

By construction, the link component L a is a periodic orbit of the local Beltrami 
field v a . Our main result in this step is the following theorem, which asserts that 
this periodic orbit is robust under suitably small perturbations. The basic idea is 
that, as our choice of the vector field w a on the strip guarantees that the cycle L a 
of the field v a has a stable manifold (which is tangent to the strip), the fact that 
the field v a preserves volume implies that the cycle also has an unstable manifold, 
thus ensuring its hyperbolicity. The proof of this result is given in Section [6j 

Theorem 3.5. Consider a tubular neighborhood N a of the link component L a , 
and an associated strip S a . Let the field v a be the unique solution to the equation 
curli> a = Xv a in this neighborhood taking the value i> a |s a = w a , where w a is the 
vector field constructed in Lemma \3.4\ Then if v is any vector field C r -close to 
the field v a in the neighborhood N a with r 1, there is a diffeomorphism $ a o/K 3 
C r -close to the identity transforming the link component L a (which is a periodic 
orbit of the field v a ) into a periodic orbit 4> Q (L a ) of the perturbed vector field v. 
Furthermore, one can assume that the diffeomorphism <£> a coincides with the identity 
outside the open set N a . 

Step 4 (Better-than-uniform approximation by global solutions). Our goal in this 
section is to present an abstract theorem that permits to approximate a local Bel- 
trami field by a global Beltrami field. The paradigmatic result on the approximation 
of local solutions of an analytic differential equation by global solutions is the cele- 
brated Lax-Malgrange theorem jT5] HH] , which applies to linear elliptic differential 
operators and yields uniform approximation on compact sets. The hypotheses of 
compactness and, most importantly, ellipticity are certainly not satisfied by the 
Beltrami equation, so our analysis requires some further elaboration. Moreover, 
better-than-uniform approximation is necessary in order to exploit the stability 
results we carefully established in Step 3. 

Better-than-uniform approximation refers to the fact that we do not only wish 
to prove that the C r norm of the difference of two vector fields is small (say, smaller 
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than a certain number e) but, rather, that the C r -approximation can be actually 
improved as one moves away from a fixed point (e.g., the origin). This is usually 
expressed using an "error function" e(x) and pointwise bounds, as we do in the 
statement of this theorem. Its proof will be provided in Section [7] 

Theorem 3.6. Let the closed set S be a locally finite union of pairwise disjoint com- 
pact subsets of K 3 whose complements are connected, and let v be a local Beltrami 
field that satisfies curlf = Xv in S. Then one can find a C r better-than-uniform 
approximation of the field v by a global Beltrami field v, i.e., for any integer r and 
any positive continuous function e(x) from S to (0, oo) there is a field v which sat- 
isfies the equation curl v = Xv in M 3 and such that its difference with the field v 
admits the pointwise C r bound 

\D a v(x) - D a v(x)\ < e(x) 

| a | ^ r 

in the set S . 

Regarding the hypotheses on the set S in the statement of the theorem, it should 
be mentioned that the condition that the complement of each compact subset be 
connected is necessary in order to apply the Lax-Malgrange theorem for an auxiliary 
elliptic equation, while the fact that the set S is the locally finite union of disjoint 
compact sets is used to attack the full problem using an iterative scheme. The 
assumption that the equation cuvlv = Xv holds in the closed set S (which means 
that it holds in an open neighborhood of S and is also a hypothesis of the Lax- 
Malgrange theorem) is required to ensure that the fields do not have a pathological 
behavior at the boundary. 

Step 5 ( Construction of the global Beltrami field). Finally, in this step we will obtain 
the desired global Beltrami field u having a set of periodic stream lines diffeomorphic 
to the link L using the local Beltrami fields v a constructed in Step 3 and the 
approximation theorem established in Step 4. By taking smaller neighborhoods if 
necessary, we will assume that the latter fields satisfy the equation curlw a = Xv a 
in the closure N a of a tubular neighborhood of the link component L a . In order 
to apply the Approximation Theorem |3.6[ we define a vector field v in the union 
S := UaGA N a of all these sets by letting it be equal to the local Beltrami field v a 
in each set iV a . Here we are using that these sets N a can be assumed disjoint for 
distinct link components, as we saw in Lemma |3.4| 

Since the field v thus defined satisfies the Beltrami equation curl?; = Xv in the 
closed set S and S verifies the hypotheses of the Approximation Theorem |3.6[ we 
infer that there is a global Beltrami field u that approximates the field v in the C r 
better-than-uniform sense. That is, given any continuous function e{x) mapping the 
set S into the positive reals (0,oo) and an integer r (the order of approximation), 
one can choose the global Beltrami field u so that in the set S we have the pointwise 
C r bound 

\D a v(x) -D a u(x)\ < e(x). 

We shall next show that the global Beltrami field u has a set of stream lines 
diffeomorphic to the link L. This will hinge on the robustness of the periodic orbits 
corresponding to each link component L a of the local Beltrami fields v a , proved 
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in Theorem |3.5| of Step 3 for suitably C-small perturbations with r 1. More 
precisely, in the present situation this theorem asserts that for any S > and any 
link component L a one can take a positive constant e a (depending on 5 and on the 
component) such that, if the vector field u is e a -close in the C r norm (with r ^ 1) 
to the local Beltrami field v a in the set N a , then there is a diffeomorphism <3> a of 
K 3 transforming the link component L a into a periodic stream line $> a (L a ) of the 
global Beltrami field u and which lies in a ^-neighborhood of the identity (that is, 
||$ a — id||c"-(R 3 ) < $)■ Since the field v coincides with the local Beltrami field v a 
in the set N a (for each labeling index a), one can certainly ensure that the global 
Beltrami field u is sufficiently close to the field v a in each neighborhood N a by using 
the aforementioned C r better-than- uniform approximation result. (More precisely, 
this amounts to taking an order of approximation r ^ 1 and an "error function" 
e(x) which is bounded from above by the previously introduced constant e a in each 
set N a .) It should be noticed that better-than-uniform approximation is essential 
for the argument, as the constants e a might not be uniformly bounded away from 
zero (with respect to the index a). 

We saw in the Stability Theorem [33] that the above diffeomorphism $ a was only 
different from the identity in the open neighborhood N a of the link component L Q , 
so these diffeomorphisms can be trivially "glued together", defining a diffeomor- 
phism $ of M. 3 as 

$ a (x) iixeN a , 

This is the diffeomorphism whose existence was claimed in the Main Theorem |l.f| 
for it transforms the link L into a set of periodic stream lines Q(L) of the global Bel- 
trami field u (which satisfies the equation curlu = Ait in R 3 ) and, by construction, 
the C r -distance between this diffeomorphism and the identity is at most 6. 

4. Proof of Theorem 13.11 

In this section we will prove Theorem |3.1| of Step 1, which ensures existence and 
uniqueness of solutions to the equation curl?; = Xv with Cauchy data v\s — w on 
the surface E (or, in terms of the associated 1-forms, *d/3 = A/3 and /3|s =7) under 
appropriate hypotheses on the vector field w. We will need to use differential forms 
of mixed degrees, so for any domain V C K 3 we use the notation 

3 

ST{V) :=0fi p (VO 

p=0 

for the graded algebra of C°° differential forms in this domain. The elements of 
$1'{V) will be denoted as 

and we henceforth identify a p-form with its natural inclusion in Q*(V). The parity 
operator Q : Q'(V) -> fi'(V) is defined by 

:= tp° © -ip 1 © ^ 2 © ~ip 3 ■ 

Clearly the differential operator curl (or its counterpart *d, which maps the space 
of 1-forms r2 1 (M 3 ) into itself) is not immediately amenable to an approach based 
on the Cauchy-Kowalewski theorem because it has a degenerate symbol. Hence, 
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we find it convenient to consider the Dirac-type operator d + d* acting on S1*(IR 3 ), 
where d* denotes the codifferential. Notice that the operator d + d* is elliptic, as 
its iterated square equals the Laplacian —A, but the price to pay for this is that, 
unlike the "curl" -kd, d + d* does not generally map 1-forms into 1-forms. 

We can invoke the Cauchy-Kowalewski theorem to derive that, for some neigh- 
borhood V of the surface S, there exists a unique C u solution ip £ il m (V) of the 
following auxiliary Cauchy problem: 

(d + d*)tp = \*ip, -01s =7, 

where 7 is the 1-form associated with the vector field w appearing in the statement 
of Theorem |3.1| as Cauchy data. 

A crucial observation is the following technical lemma, which asserts that the 
hypotheses made on the Cauchy data imply that the unique solution of the above 
problem is co-closed (i.e., d*tp — 0). This will be used to relate solutions of the 
above auxiliary Cauchy problem to solutions of the equation curlu = Xv, v\s = w: 



Lemma 4.1. Under the hypotheses on the 1-form 7 imposed in Theorem \3.1\ 
(namely, that its pullback to the surface is closed and that its dual vector field 
is tangent to the surface), the solution of the Cauchy problem (d + d*)ip = A*ip 
with ip\s = 7 is co-closed, i.e., d*ip = 0. 

Proof. By the identity *<i* = —Qd*, taking the action of *d on the equation (d + 
d*)ip — A * ip we readily obtain that the form d*ip satisfies the PDE 

{d + d*)d*ijj = -\*Qd*ip. 

The Cauchy-Kowalewski theorem will then yield that the codifferential of the form 
ip is identically zero in its domain V provided that it vanishes on the surface S. 

In order to show that the codifferential d*ip is zero on £, let us use the local 
coordinates (p, £ x , £ 2 ) introduced in Step 1 of Section|3]and denote by J the exterior 
derivative in the variables (^ 1 ,^ 2 ), so that the hypothesis that the 1-form ^(7) is 
closed is equivalent to the condition djl^: = 0. Let us decompose the form ip as 

ip =: * + dp A $ 

for some forms ^f,^f £ fl*(V) without dp (by which we mean that the interior 
product of the coordinate vector field 8/ dp with them vanishes identically). As 
"01s — 7 is a 1-form without dp and one has d'yls — by hypothesis, it is easy to 
check that, on the surface, the star of the form ip can be written as 

= *7|s =: dp A 7 
for some 1-form 7, while its differential reads as 

il , 5* 

ii\V 1 = Hn A 

dp 

The codifferential of the form ip on the surface can then be computed using the 
equation (d + d*)ip — A * ip, finding that 

(4.1) dV|s = (A*^-#)| E = dpA (A7- . 

The degree-0 component of the right hand side of the above equation is necessar- 
ily 0, which implies that the codifferential of the degree-1 component of ip vanishes 



dip\i 



d^> + dp A I — d^> 
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on the surface, i.e., cZ* Is = 0. For degrees p = 2,3, the fact that the form ip p is 
zero on the surface S implies that 



d*(r)U = (-1) P (*M P )| S = (-17* (dp A 



d*ip p 



dp 

which is a (p— l)-form without dp. Comparison with the degree-(p — 1) component 



of d*tp, computed in Eq. (4.1), then shows that d*(V> p )|s = for p = 2,3 and the 



lemma follows. □ 
As the form of mixed degree r/> solving (d+d*)tp — X-kip with t/>|a = 7 is co-closed 



by Lemma 4.1 above, from the identity *(*^) = ip we deduce that the form ip also 
satisfies the equation -kdtp = Xip- Since -kd maps 1-forms into 1-forms, the degree-1 
component j3 :— ip 1 is then the unique solution of the Cauchy problem 

*dp = \p, /3|s=7 



in a neighborhood of the surface, and hence Theorem 3.1 follows by letting v be 



the vector field associated to the 1-form j3 via the Euclidean metric. 



5. Proof of Lemma l3~4l 

In this section we construct a vector field w a , tangent to a strip E a contained in 
a tubular neighborhood N a of the link component L a , which satisfies the hypothe- 
ses of the Existence Theorem |3.1| and whose pullback to the surface has the link 
component as a hyperbolic limit cycle. As we will see, the 1-form j a dual to the 
desired field w a is simply given by dO — z dz in a natural coordinate system; hence, 
we will start by carefully defining these coordinates on the neighborhood N a . The 
coordinates we shall use are (p, z, 9), where p is the signed distance to the surface 
E a and, roughly speaking, z and 9 are suitable extensions of the signed distance to 
the link component L a as measured on the surface and of the arc-length parameter 
of L , respectively. 

To construct the coordinates, we begin with an arc-length parametrization of the 
link component L a , which is a periodic function 9 mapping the component to the 
reals modulo the length \L a \ of the link component (i.e., R/|L n |Z). By taking the 
strip small enough, we can safely assume that the signed distance function to the 
link component L a , as measured along the surface E a with respect to the metric 
h inherited from its embedding in M 3 , is analytic. Let us denote this function by 
z : E a —> M. The gradient operator associated to the metric h on the strip will be 
denoted by V. 

We can extend the parametrization 9 of the link component to a C u function 
on the surface, which we still call 9 with a slight abuse of notation. Calling (f> t the 
local flow of the gradient field Vz in the surface E a , this extension is defined by 
setting 

9{4> t x) := 9{x) 

for all points x in the link component L a and all t for which <p t x is a well-defined 
point in the surface E a . By periodicity, this extension takes values in the reals 
modulo the length \L a \. Obviously V9 7^ in the surface E a , and for any point 
y G E a there is a unique pair (x, t) e L a x K such that y = <fi t x (in fact, here t is the 
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signed distance z(y) because h(Vz, Vz) — 1). Since the gradients of the functions 
9 and z are orthogonal by construction, i.e., 

(5.1) h(V6, Vz) = , 

it is apparent that [z, 9) is a coordinate system on the surface E . 

To construct the adapted coordinate system in the tubular neighborhood N a , 
let us begin by considering the signed distance function to the surface E a in N a , 
which we denote by p. This function is analytic if the neighborhood N a is small 
enough [T3]. We now extend the coordinate z (resp. 9) on the surface to a C u 
function mapping the neighborhood N a into the reals (resp., by periodicity, the 
reals modulo the length \L a \); for simplicity, we still denote the extended functions 
by z and 9. This extension is performed using the local flow <f> t of the gradient field 
Vp as 

z(4> t x) := z(x), 9{<j> t x) := 9(x) , 
for all points x £ S Q and values t such that <fitx £ N a is defined. Arguing as above, 
it can be readily established that (p,z,6) is an analytic coordinate system in the 
open set N a verifying the conditions 

(5.2) (Vp,V9) = (Vp,Vz)=0, \Vp\=l. 

It should be noticed that the surface S a is given by the zero set {p = 0} and the 
link component L a is {p = z = 0}, which implies the relations 

(5.3) (Vz,W)L =0, |V0|L =1 



as a consequence of the orthogonality condition (5.1) and the fact that 9\^ a is an 
arc-length parametrization of the link component L a . 

As we mentioned, in this coordinate system we can define the 1-form associated 
to the desired field w a as 

"f a '■= d9 — z dz . 

The pullback to the strip of this 1-form is obviously closed. The associated vector 
field, w a := V# — z Vz, is clearly tangent to the strip E Q because it is orthogonal 
to Vp by Eq. (5.2). Moreover, by the relations (5.3) it follows that the field w a is 
nonvanishing and orthogonal to Vz on the link component L ai which implies that 
this component is a periodic orbit of the field w a . 

That the link component L a is an asymptotically stable periodic orbit of the 
pullback ja(w a ) of the field w a to the surface will stem from the existence of an 
appropriate Lyapunov function. This Lyapunov function is simply z 2 , since it is 
nonnegative, vanishes exactly on the cycle and its scalar product with the field w a , 



( Wa ,V(z 2 ))| Sa =2(z(Vz,V0)-z 2 |Vz| 2 )L =~2z 2 |Vz| 2 



is zero on the cycle L a — {z = 0} D S a and strictly negative in its comple- 
ment ^ a \La- 

To prove the hyperbolicity of the cycle L a , we start by observing that the field 
w a on the surface has the form 

f « {Wa) = g{z > e) t-9~ Z iz 



on account of the relations (5.2) and (5.3), where the function g{z,9) is defined in 
terms of the gradient of 9 and the metric h on the surface as g := h(\79, V0). Hence, 
the integral curve T(t) of the field ja{w a ) corresponding to the link component L a 



KNOTS AND LINKS IN STEADY SOLUTIONS OF THE EULER EQUATION 



13 



is parametrized by z(t) = and 9{t) = t modulo the length \L a \. The monodromy 
matrix along this orbit has only one nontrivial eigenvalue p, which by Liouville's 
formula (cf. e.g. [5J Theorem IV. 1.2]) can be expressed as 



\x = exp 



= e -|£.| <i 



using that the function g o T is |L a |-periodic. Hyperbolicity then follows. 

6. Proof of Theorem 13.51 



We shall prove the Stability Theorem 3.5 which ensures that the link compo- 



nent L a is a robust periodic stream line of the local field v a , which solves the 
equation curlu a = Xv a with Cauchy data v a \^ a = w a chosen as in Step 2. To do 
this, we will use the same coordinate system (p, z, 9) we introduced in the proof of 



Lemma 3.4 defined in a tubular neighborhood N a of the cycle. 



From the construction of the vector field w a used as Cauchy data (cf. Lemma 3.4 ), 
it stems that the local Beltrami field v a can be written as 

v a = (g(9) +0{z)+ 0{p)) ^-(z f(9) +0(p) + 0(z 2 )) ± + 0{p)-^. 



In fact, from the expression for the pulled back field ia{w a ) given in Eq. (5.4) it 
follows that we can set f{9) = 1 above. As the field v a is analytic, the small terms 
behave like symbols in the sense that, e.g., 



jO{p s ) 



< F mnl {z,9)\p\ s 



dp n dz m d9 l 

for some smooth function F mn i. The function g(6) entering the above expression 
for the field v a is nonvanishing and |L a |-periodic, \L a \ being the length of the link 
component, which means that we can rescale the coordinate 9 to have g(9) = 1. 
The integral curve T(t) of the field v a corresponding to the link component L a , 
which will have certain period T after the rescaling, is given in these coordinates 
by p{t) = 0, z{t) = and 9{t) = t modulo T. 

A key intermediate result is that the cycle L a is hyperbolic. The proof is based 
on a monodromy matrix computation: 

Lemma 6.1. The link component L a is a hyperbolic periodic orbit of the local 
Beltrami field v a . 

Proof. The variational equation associated with the integral curve T(t) correspond- 
ing to the link component can be written as 

f t =Dv a (T(t))y(t), 

where y(t) is a three-component vector and 

D««(T(t)) = F 2 (t) 
W*j 

is the Jacobian matrix of the field v a in the coordinates (p, z, 9) evaluated at the 
integral curve T(t). Here Fj(t) (1 ^ j ^ 4) are some analytic T-periodic functions 
whose explicit expression we will not need. 
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The solution of the variational equation above is given by y(t) = M(t)y(0), 
where 

/Gi(i) 0\ 
M{t) := G 2 (t) e-* 
\G 3 (t) G 4 (t) 1/ 

is the fundamental matrix and Gj(t) are analytic functions. By definition, the 
monodromy matrix is M(T), with T being the period, so its nontrivial eigenvalues 
are \x\ := G\(T) and /l«2 := e _T . 

The field v a is divergence- free because it satisfies the equation curliv = \v a , 
so a direct application of the Liouville formula shows that the determinant of the 
monodromy matrix is given by 

/ii^ 2 = exp f J (div v a ) o T(t) dt \ = 1 . 

As the eigenvalue fi 2 = e _T is smaller than 1, we deduce that the other eigenvalue 
/ii is bigger than 1, thus establishing the hyperbolicity of the cycle. □ 

Let us now conclude the proof of the Stability Theorem |3.5| Since the link compo- 
nent L a is a hyperbolic periodic orbit of the field v a , a straightforward application 
of the hyperbolic permanence theorem [TTJ Theorem 4.1] shows that any vector 
field v sufficiently C r -close to the field v a in the neighborhood N a has a periodic 
orbit which is a C r -small deformation of the component L a , provided r ^ 1. More 
precisely, for any integer r ^ 1 and any S > there exists some e > such that any 
vector field v with \\v a — v\\q, -r N \ < e has a periodic orbit A in the neighborhood 
N a , and moreover the cycles L a and A admit smooth trivializing maps O a and 
from N a to M 2 such that the orbit L a (resp. A) is the zero fiber of the map a (resp. 
of 6) and the difference of the latter maps is C r -bounded as ||0 O — Q\\c r (N a ) < 
The existence of the desired diffeomorphism $ a of M 3 then follows from Thorn's 
isotopy theorem [H Theorem 14.1.1]. 

7. Proof of Theorem 13.61 

In this section we shall prove Theorem |3.6| which guarantees that any local 
Beltrami field v defined in a closed set S satisfying certain technical conditions can 
be approximated in this set by a global Beltrami field v in the C better-than- 
uniform sense. Let us denote by Sb the connected components of the closed set S, 
with b ranging over an at most countable set B. Since the set S is the locally finite 
union of Sb, it is clear that there exists an exhaustion =: Kq C K\ C K2 C • • ■ 
by compact sets of M. 3 such that: 

(i) The union of the interiors of the sets K n , with n € N, is M 3 . 

(ii) For each n, the complements of the sets K n and of S U K n are connected. 

(iii) If the set K n meets a component Sb of S, then Sb is contained in the 
interior of K n+ i. 

A key technical lemma is the following better-than-uniform approximation result 
for the auxiliary scalar elliptic equation (A + A 2 )/ = 0, which is modeled upon a 
theorem of Bagby and Gauthier [2] . The proof is based on an indirect induction ar- 
gument that uses the Lax-Malgrange theorem and the properties of the exhaustion 
we have just introduced. 
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Lemma 7.1. Let the closed set S be a locally finite union of pairwise disjoint 
compact subsets of M 3 whose complements are connected. If the scalar function f 
verifies the equation (A + A 2 )/ = in the set S, then one can find a C s better-than- 
uniform approximation of the function f by a global solution g of the aforementioned 
equation. That is, for any integer s and any positive continuous function e(x) from 
the set S to (0, oo) there is a function g satisfying the equation (A + \ 2 )g = in 
M 3 and such that its difference with the function f has the pointwise C s -bound 

\D a f(x)-D a g(x)\<e(x) 

in the set S . 



Proof. The proof relies on an induction argument that is conveniently presented in 
terms of a sequence of positive numbers {e n )^ =1 related to the error function e{x) 
and to the above exhaustion by compact sets K n through the conditions 



(7.1) 



1 

< — mm 
6<7 xeK m4 



e(x) and 



oo 

E ■ 

n—m+1 



< e r , 



which are required to hold for all m 1. For later convenience, we set eo : — and 
we are dividing by the number a := card{a : \a\ ^ s} of multiindices a with \a\ 
smaller than or equal to the order of approximation s. 

The induction hypothesis is that there are functions g n : M 3 — > K satisfying the 



equation (A + \ 2 )g n = in 
C s estimates hold: 



and such that, for any integer p 1, the following 



(7.2a) 

(7.2b) 
(7.2c) 



sup 

Sn(K p+1 \K p ) 



sup 

Sn(K p \K p - 



D a 



D a 



p 

n=l 
P 



=1 



< e r . 



2c 



sup 



\D a g P \ < e p 



p-i 



£p-l 



Here and throughout the proof of this lemma, all the multiindices will be assumed 
to range over the set \a\ ^ s without further mention. 

Let us start by noticing that, by the Lax-Malgrange theorem jTS], [T5], there 
exists a function 171 : M 3 — > K satisfying the equation (A + X 2 )gi = in M 3 and 
the C s estimate \D a (f — g\)\ < t\ in the set S H K%- Since we chose Ko = and 
eo = 0, it is a trivial matter that the induction hypotheses (7.2) hold for p = 1. We 
shall hence assume that the induction hypotheses hold for all 1 ^ p ^ m and use 
this assumption to prove that they also hold for p = m + 1 . 

To this end, let us construct a function f m on the set SU K m by setting f m := 
in the set K m and defining f m on each component Sb of the set S as 



/m| St, '■— 



f ~ J2 9n 

n=l 





if s b n (K m+2 \ K r. 
if S b n (K m+2 \ K r. 



0* 



-1) = 



Here K n stands for the interior of the set K n . The definition of the exhaustion 
and the first induction hypothesis (7.2a) guarantee that the function f m satisfies 
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the equation (A + A 2 )/ m = in its domain and that one has the C s estimate 



n=l 



(7.3) sup \D a f m \ sC sup 

K m U(SnK m + 1 ) sn(K m+1 \K m ) 

A further application of the Lax-Malgrange theorem allows us to take a function 
.9m+i which satisfies the equation (A + X 2 )g m+ i = in R 3 and which is close to the 
above function f m in the sense that 



(7.4) 



sup 

K m+2 n(SUK m ) 



\D a {fm — 9m+l 

)|< 



Eq. ( 7.4 1 above and the way we have defined the function f m in the set K m+ 2\K m+ i 
ensure that the first induction hypothesis ( 7.2a I also holds for p = m+l. Moreover, 
from the relations (7.2a), ( |7.3| and (7.4 1 one finds the following pointwise estimate 
in the set S D (K m+ i\K m ): 



m+l 



D a [f 



71=1 



\D a g, 



m+l 



< e. 



-i)| + |D a /J < £m+i + 2e r 



This proves the second induction hypothesis ( |7.2b ) for p = m + 1. Furthermore, 
sup |-D a g m+ i| < sup \D a (f m — g m +i)\ + sup |D a / TO | < 



by the relations (7.3) and (7.4), so the third induction hypothesis (7.2c) also holds 
for p = m + l. This completes the induction argument. 
The desired global solution g can now be defined as 



oo 

E 

n=l 



fin 



with this sum converging C s -uniformly by the definition of the constants e n (see 
conditions (7.1)) and the third induction hypothesis (7.2c). As the functions g n 
verify the equation, it is easily checked that the function g also satisfies the equation 
(A + X 2 )g = in R 3 , which in turn ensures that g is analytic by elliptic regularity. 
In addition to this, from the definition of the constants (7.1) and the induction 
hypotheses (7.2) it follows that, for any integer m, in the set S D (K m+ i\K m ) we 

-\D a g m+2 \+ D c 

< cr [(e m+ i + 2e m ) + (e m+2 + e m +i) + 



have the pointwise C s estimate 

D a 



\D a (f-g)\ < 

|a|<s 



E 



m+l 

/ ~ E 9n 
n=l 



E 9n 

n— m+3 



E 

n— m+3 



(e„ + e n -i)] 



< 2a(e m + 2e TO+ i) < min e(x) , 

xeSn(K m+1 \K m ) 



where a was the number of multiindices a with |a| ^ s. The better-than- uniform 
approximation lemma then follows. □ 



Let us now complete the proof of the Approximation Theorem |3.6[ As the 
local Beltrami field v satisfies divi> = in the set S, acting with the operator 
curl+A on the equation curlw = Xv we readily derive that (A + X 2 )v = 0, i.e., 
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that each component Uj of the vector field (in Cartesian coordinates) satisfies the 



scalar equation (A + A = 0. By Lemma 7.1 above, we can then approximate 
the function Vi by a global solution wi of the equation (A + X 2 )wi = in the C s 
better-than-uniform sense. To be more concrete, for later convenience we shall take 



s = r + 2 (with r the order of approximation in Theorem 3.6) and choose an error 
function e(x) to be specified later, thereby assuming that the global solution Wi 
approximates the function Vi pointwise in the set S as 



E 

|a|<r+2 



\D a Vi(x) - D a Wi(x)\ <e(x), xeS. 

Let us now denote by w the vector field in E 3 whose Cartesian components are 
Wi. Defining the auxiliary vector field w :— curl wj A, an easy computation using 
the above C r+2 estimate for the difference Vi — Wi yields that the difference of the 
fields v and w satisfies the following C r+1 pointwise bound in the set S: 
(7.5) 

\D a v(x) - D a w(x)\ = i \D° curlv(x) - D a cur\w(x)\ < ^1 . 

\a\^r+l \a\^r+l 

By construction, the vector field w obviously satisfies the equations (A + X 2 )w = 
and div w = 0, which in turn implies that 

(curl +A) (curl- A) w = 

because, being divergence-free, the Laplacian of this field is Aw = — curl curl w. 
Hence the analytic vector field v := (curl +A)w/(2A) satisfies the Beltrami equation 



curlw = \v and, as a consequence of the estimate (7.5), the pointwise C r bound 



£ \D<*v(x)-D"v(x)\ = ± Yl \D a (cuA+X)(w(x)-v(x))\< (6 + 

in the set S. Thus v is the field whose existence was claimed in the theorem, 
provided that the error function e(x) is chosen small enough (e.g., smaller than 



A e{x) /(6 + 3A), where e(x) is the error function in the statement of Theorem 3.6 ) 



Remark 7.2. It should be noticed that we have used that the parameter A is nonzero 
in the last part of the proof; in fact, if A — one cannot hope to approximate the 
local solutions by global ones because an irrotational vector field defined in a simply 
connected region cannot have a periodic trajectory. 



8. Applications and examples 

To illustrate the Main Theorem we shall address a problem of Etnyre and 
Christ concerning the existence of steady solutions of the Euler equation having 
stream lines of all knot and link types. We shall also indicate how Theorem |1.1| can 
be extended to any open analytic Riemannian 3-manifold. 

Example 8.1. While in the literature one can explicitly find questions regarding 
the existence of solutions whose stream or vortex lines are simple links (e.g., the 
Borromean rings in [2D] ) , a particularly attractive question on the complexity of 
steady Euler flows was formulated by Etnyre and Ghrist in [6] : Are there any steady 
solutions to the Euler equation having stream lines of all knot and link types in 
Euclidean 3-space? 
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Figure 1. Example of a locally finite link containing the Bor- 
romean rings, the trefoil, the figure eight and the (7, 4) knots. 

Theorem 1 1 . 1 1 readily provides a positive answer to Etnyre and Christ's question 
for locally finite links. To see this, it suffices to notice that there is a countable 
number of isotopy classes of tame knots and locally finite links in K 3 [3T], which 
allows us to construct a locally finite link L containing a sublink of each isotopy 
class by the noncompactness of M 3 (cf. Figure [l] for an example of this kind of 
construction). 

Example 8.2. A straightforward modification of the proof of Theorem [O] actu- 
ally works for the Euler equation in a general Riemannian open 3-manifold M of 
class C w . Let us briefly elaborate on this point. The Euler equation automatically 
makes sense on a Riemannian manifold, with V being the covariant derivative as- 
sociated to the metric g and u a vector field in M. A Beltrami field is most easily 
characterized in this context in terms of the associated 1-form j3, which satisfies the 
equation -kd/3 = A/3. 

Mutatis mutandis, both the statements and proofs in the Steps 1 to 3 of the proof 
of Theorem 1 1 . 1 1 remain valid in this more general situation. This is not surprising, 
as they correspond to a (rather delicate) local construction which does not make 
use of any particular feature of the Euclidean metric. 

Step 4 is also valid after minor modifications. In the non-Euclidean setting, it 
is convenient to work directly with the 1-form /3, so that the role of the Euclidean 
Laplacian in Section [7] is now played by the Hodge Laplacian A := —(dd* +d*d). It 
is evident that better-than-uniform approximation of vector fields can be painlessly 
translated into better-than-uniform approximation of the associated 1-forms (gen- 
erally with different error functions) . The 1-form j3 then satisfies the equation 



(A + X 2 )f3 = 0, for which Lemma 7.1 and its proof can be easily adapted. The 
remaining parts of Step 4 remain unchanged, and so does Step 5. The analog of 
Theorem II .11 on M thus follows. 
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